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I. INTRODUCTION 

A full adder is a basic component of a conventional computer and a welcome asset for 
quantum computers. In particular, the Shor quantum algorithm requires modular exponen- 
tiation, f[x) = a'^(modA^), which cannot be computed without a quantum full adder. The 
question arises: What is a full adder in a quantum computer? A quantum computer oper- 
ates on a superposition of numbers simultaneously. What is not clear is what it means to 
add two superpositions. One definition of a full adder in a quantum computer is that a full 
adder is a gate which adds a given number to a superposition of numbers. This full adder 
must simultaneously add a definite number, "a" , to all numbers, , which are in a super- 
position in a quantum computer register. The basic idea for the full adder is well-known 
However it is not clear how to implement this idea in a many-qubit spin quantum 
computer. It is even less clear how to simulate accurately a quantum full adder having a 
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large number of qubits using a conventional computer. One must understand the role of 
non-resonant effects and how to minimize them. One must also know the structure of the er- 
ror states caused by non-resonant effects. Our paper provides an answer for these issues. To 
implement a quantum full adder, we propose to use a sequence of electromagnetic vr-pulses 
on a spin-chain quantum computer. We have simulated a dynamics of this quantum full 
adder with 201 spins, on a conventional computer; analyzed unwanted non-resonant effects; 
and determined the structure of the error states and ways to reduce non-resonant effects. 



A classical full adder operates with an input of two addend bits, "a" and "b", and a 
carry bit, "c". (See Table 1.) In Table 1, s' and c' are the output sum and the carry-over, 
respectively. The sum, s', can be easily expressed as s' = a © 6 © c (where © is an addition 
modulo 2). The formula for a carry-over is more complicated. One can see from the Table 
1, that c' = ab (B ac(B be. 



We describe the basic idea of a quantum full adder first suggested in [Q. The full adder 
quantum gate (F) depends on the value of the bit "a". If a = 0, the quantum computer 
applies the gate F{0). If a = 1, it applies F{1). The F-gate can be expressed in terms of 
the Control-Not (Cik) and Control- Control-Not (Cikp) gates. A Cj^-gate changes the value 
of a target qubit, "k", if the control qubit, "i", has the value "1". A Cj^p-gate changes the 
state of the target qubit, "p", if both qubits, "i" and "k", have the value "1", 
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Table 1: Table of binary addition. 

We shall also use the Not- gate which can be designated by It changes the value of the 
ith qubit independent of the values of all other qubits. Using these gates, we can transform 

the state, IO3C261), into the state, IC3S261), 

F32i(a) ^ IO3C261) = 14461), (2a) 

where "a" and "b" are the addend bits and s' is their sum, c and c' are the input and output 
carry bits. Namely, 

-^321(0) = C12C123, -^321(1) = C'l2C'l23C'2C'23- (26) 

We use the convention that the right gate acts first. Let us check the action of the full adder, 
F321, for example, for a = 6 = c = 1. We use ^321(1): 

C23IO3C261) = C23IO3I2I1) = II3I2I1). (3) 

In (3), the second control qubit has the value "1". So, the third target qubit changes its 
value from "0" to "1". Next, 

C2II3I2I1) = |l302ll),Ci23|l302ll) = II3O2I1), (4) 
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C12I130211) = I131211) = \cw2bl). 

Thus, in accordance with the Table 1, we obtain, 

b=l, s' = 1, c = 1. 

IV. A SPIN CHAIN QUANTUM COMPUTER 

An atomic chain quantum computer based on triplets, ABC ABC ABC..., was first in- 
troduced in P). The implementation of this idea for a chain of spins interacting through the 
Ising interaction was given in 0]. Ising-type spin chains have been used for quantum com- 
putation in a statistical ensemble quantum computers [^|[. In [0, we considered a simple 
model - a chain of identical spins in a non-uniform magnetic field. 

The present paper is based on and [0. We consider a chain of spins (e.g. nuclear spins) 
in a non-uniform magnetic field. Similar to |]^, the angle, B, between the direction of the 
chain and the direction of the permanent magnetic field (5;-direction) satisfies the condition: 
cos = 1/ \/3. Then, the dipole-dipole interaction between spins is suppressed, and the Ising 
interaction becomes dominant. We suppose that our chain consists of a periodic sequence 
of triplets, ABC ABC ABC .... The triplet, ABC, can be different nuclear spins or identical 
spins in slightly different environment. For definiteness, we shall keep in mind this second 
case with the following typical parameters, 

cji/27r ~ 400Mifz, Uk = uji + {k - l)Au;, Aa;/27r ^ 2{]kHz, (5) 

Jac/'^t^ ~ lOOifz, Jbc = '^Jac, Jab = ^Jac, 

where Ui is the Larmor (angular) frequency of the right end spin, Uk is the Larmor frequency 
of the kth spin, Jj^ is the constant of Ising interaction between the ith and kth spins. 

In the presence of a circularly polarized transverse magnetic field, the Hamiltonian of 
the system can be written as 



k 

where P is the nuclear spin operator, cu and (f are the frequency and the phase of the 
transverse magnetic field, ft is the Rabi frequency (the magnitude of the transverse magnetic 
field in frequency units), h = 1. 

V. IMPLEMENTATION OF THE QUANTUM FULL ADDER 

First we shall define our problem. Suppose that we have a number \b) — \b^^~^\..b^°^), 
J)M = 0, 1. In decimal notation, b = b^^^^^2^^^ + ... + 6*^°-' 2*^. (Below, we omit parentheses in 
the superscripts.) Note, that in general, we have a superposition of many numbers, \b)i, and 
any gate must act on all of these numbers simultaneously. We are going to add to a number 
\b), (all numbers, bi, in a superposition) a definite number, a = a^^^...a^, where a™ = 0, 1. 
To achieve this goal we shall use 2L + 1 qubits. We load the number 6 in a chain of qubits 
in the following way, 

\02L+ib^E'02L-ib^L-2-05b%02b',). (7) 

This means that we place two additional qubits in the states |0) in front of the qubit b^, and 
one additional qubit in the state |0) in front of all other qubits fe'", m 7^ 0. 

The gate C123 is not convenient for our spin chain quantum computer in which each spin 
interacts only with its neighbors. So, we replace it with the following sequence of gates, 

C123 = C*23C'32C'23C'l32C'23C'32C23. (8) 

Let us check this equation, for example, for the state, IO3I2I1). Prom the left side of Eq. (8) 
we have, 

C123IO3I2I1) = II3I2I1), (9) 

as the first and the second qubits are in the states, the third qubit changes its state. 
Next, we follow the operations on the right hand side of (8) to show that the same result is 
obtained. Prom the right side of Eq. (8) we have, 

C23IO3I2I1) = II3I2I1), C32II3I2I1) = II3O2I1), (10) 
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C23II3O2I1) = |l302ll>, Ci32|l302ll> = II3I2I1), 
C23II3I2I1) = IO3I2I1), C32IO3I2I1) = IO3I2I1), 
C23IO3I2I1) = II3I2I1), 

which coincides with the right side of Eq. (9). Now, instead of (2b), we have the following 
expression for the full adder, F, 

-^321(0) = C12C23C32C23C132C23C32C23, (11) 

-^321(1) = C'l2C'23C'32C'23C'l32C*23C*32C23C2C'23- 

Now we explain how to add the numbers b and a. If a° = 0, the quantum computer applies 
the gate, -^321(0). If a° = 1, the quantum computer applies the gate, ^321(1). According to 
(2), the result is, 

F321IO3O26;) = 145^6?), (12) 

where we have replaced s' by s° = a° © 6°, and c' by which is the carry-over for the next 
addition of and a^. 

Next, consider the five right most qubits, 

\0,blclsy,). (13) 

To add b^ and and the carry-over, c^, we should first make a swap, S, of the values of the 
fourth and the third qubits, 

S,s\blcl)^\clbl). (14) 
The swap gate, Sik, can be represented in terms of Cik gates, 

Sik — CkiCikCki- (15) 

Let us check, for example, the action of the 5'2i gate on the state, IO2I1). Using (15), we 
have, 

'S'21 = C12C21C12, (16) 
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C12IO2I1) = |l2ll), C21II2I1) = II2O1), C12II2O1) = II2O1). 

Thus, the ^-gate transforms the state, |01), into the state, |10). 
After the action of the 5'43 gate, the state (13) changes into, 

MlsX)- (17) 
The state, IO5C463), has the form (2a), and it is ready for apphcation of the full adder, F543, 

F543\05c\bl) = \cls\bl). (18) 

Certainly, we use ^543(0) if = and ^543(1) if — 1. Thus, we obtain the sum, 
— ©fe-*^ ©c^, and a carry-over, c^, for the next addition. It is clear that by repeating the 
application of gates F and S, we shall obtain the desired answer. A complete full adder, F, 
can be represented as a combination of the elementary full adders, F^k, and the swap-gates, 
Sij, 

F = F2L+l,2L,2L-l{a^ ^)S2L,2L-1--- (19) 
i^765 (a')565i^543 (0^)^43^321 (a°) . 

After the action of the full adder, the superposition of the states (7) transforms into the 
superposition of states, 

|c^s^-i6^-is^-26^-l..s060^. (20) 

Thus, the qubits in even positions and the left end qubit carry the results of addition. The 
problem is: How to implement all of these gates using 7r-pulses? We suppose that fl ^ Jik, 
say Q/27r fa lOHz. This allows us to excite each spin individually. To implement F321, we 
use equations (11). As an example, we show how to implement ^321(0). According to Eqs 
(11), first we implement C23, then C32, and so on. Suppose we have an integer number of 
triples, ABC, in our spin chain. Then, we have 4 possible frequencies in the position (3n — 2) 
which correspond to a spin C: 

'^3n-2 — ^3n-2 + J BC + JaC, (21) 
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^3n-2 — ^3n-2 + JbC — J AC, 
'^3n-2 — ^3n-2 — J EC + J AC, 
'^3n-2 = '^3n-2 — J EC — JaC- 

Here corresponds to the kth spin whose left neighbor is in state \i) and whose the right 
neighbor is in state \ Similar expressions can be found for spins in positions (3n — 1) and 
(3n), the B and A spins, respectively. The end spins have only two frequencies. For the left 
end spin, A, 

"^^L+l = '^2L+1 + JaE, ^2L+1 = <^2L+1 — JaB, (22) 

and for the right end spin, C, 

ul^ui^ Jbc, Jbc- (23) 

Now, to implement C23, we apply a 7r-pulse with frequency uj*^ and then a 7r-pulse with 
the frequency uj\^. One of these two pulses definitely changes the state of the third spin if 
the second spin is in its excited state, To implement the gate C123, we need a single 
TT-pulse with frequency uj^. To implement a Not-gate, C2, which appears in ^132(1), we 
have to apply four tt- pulses with all possible frequencies, 0J2, i = 0, 1, and j = 0, 1. The 
total number of pulses required to implement a Fjjfe(0)-gate (if i 7^ 21/ + 1, A; 7^ 1) is 15: two 
pulses for each Control-Not gate and one pulse for the Control-Control-Not gate. For the 
Fijk{l) gate, the total number of pulses is 21 (16 pulses for the left triple of the chain). The 
swap gate, Sij, requires 6 pulses. 

Thus, to add L-qubit numbers, our proposal requires (2L -|- 1) qubits and less than 27L 
TT-pulses. 

VI. SIMULATIONS OF A QUANTUM FULL ADDER 

For numerical simulations of a quantum full adder we used the following assumption: the 
frequency difference between two neighboring spins is much greater than the Rabi frequency. 
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As a result, the selective excitation of a chosen spin has a small effect on all other spins. 
At the same time, we take into account the action of a vr-pulse on non-resonant states. For 
example, suppose that the frequency of a vr-pulse is: uj = ul^. This pulse is resonant with 
spin "3" only if the neighboring spins, "1" and "2", are in their excited states. The states 
in which the spin "3" has frequencies, u^^, u^^, and u^^, are non-resonant for the pulse with 
u = ul^. We take into consideration the transformations of all these non-resonant states. 

Note, that any pulse in our simulations acts on all basic states in the quantum superpo- 
sition. The resonant state transforms into a state with the opposite direction of the resonant 
spin. Every non-resonant state transforms into the superposition of two states: the initial 
one and an error state generated by the non-resonant transition. If the probability of an 
error state is less than a chosen small number, e, our computer program automatically re- 
moves it. One can argue that every removed state can generate a number of new states. 
However, suppose that we have removed a state with a small probability, P. Due to the main 
property of the unitary transformations, the sum of the probabilities of all states generated 
by the removed state (including itself) is P. Thus, the total probability of all neglected 
states cannot increase in spite of generation of new states! 

We have simulated the addition of the 100-digit numbers in a 201 spin chain. To minimize 
non- resonant errors, we used the 27r/i;-method |0J^. The basic idea of the 2-7? /c- method is the 
following. One chooses the Rabi frequency, Q, of a resonant vr-pulse in such a way that it 
becomes approximately a 2vrA;-pulse for all non-resonant transitions (where k is an integer 
which generally is different for different states). A 2vrA;-pulse does not generate unwanted 
error states. A transformation of a basic state under the action of a vr-pulse was described in 
our previous paper 0, where the following values of dimensionless parameters were chosen, 

Jac = I5 Jbc = 2, Jab = 3. (24) 

We found that all non-resonant transitions approximately satisfied the 2vr/c-conditions for 
the value of the Rabi frequency, Q = Qq = 0.10005. 
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FIG. 1. The probability of the expected state, Pq; and the number of error states, A^, as a 
function of the Rabi frequency, ft. 

Next, we present the results of our numerical simulations. As an example, we add the 
"classical number" , 

l9Vl..0°, (25) 
(2^^ in a decimal notation) to the "quantum number". 



(26) 



(1 in a decimal notation). The sum of these two numbers. 



(27) 



corresponds to the quantum state of the chain of 201 spins. 
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FIG. 2. "Line spectrum" of the probabilities, P, of error states (fi = 0.10021); (a) the region: 
P ~ 10"^, (b) the region: P ~ 10-^°. 

For the value of fl — flo which approximately satisfies 27r/c-condition, we have for the 
probability of the expected state (28): Pq — 0.99889. The number, N, of error states with 
probability P > 10"^^ is only 304. Small deviations from the 27rA;-condition significantly 
infiuence the result. As an example, for O, — 0.10021, the value of Pq decreases to 0.98300, 
and the number of error states, N, grows to 46530. Fig. 1 shows the dependence of Pq 
and on the Rabi frequency, Q in the vicinity of Qq. One can see that the probability, 
Pq, smoothly decreases by approximately 1% when the Rabi frequency, fl, shifts from flo 
by approximately 0.1%. Unlike the value of Pq, the number of error states with probability 
> 10~^^ does not change in the close vicinity of Qq. But it sharply increases when deviation 
of fl from flo approaches 0.15%. Fig. 2 shows the probabilities of error states (the states 
are shown in the order of their generation). One can see a specific "line spectrum" of the 
probabilities: the error states arc "attracted" to a few discrete values of the probability. 
Very similar effects were obtained for other examples of addition including the addition of 
a "classical" number to a superposition of two quantum numbers. 



11 



VII. PHASE CONTROL 



The quantum full adder implementation considered above provides a transformation of an 
arbitrary superposition of "quantum numbers", g^, into a superposition of numbers (qk + d), 
where a is any given "classical number". However, this proposed scheme also generates 
complicated phase differences between the states of the quantum superposition. This effect 
can be inappropriate, especially for the Shor algorithm. In this section, we describe how to 
extend our simple model to incorporate phase restoration after the action of every 7r-pulse. 

Consider the chain of paramagnetic ions ABCABC containing nuclear spins in a high 
nonuniform magnetic field. Every electron spin interacts with the nuclear spin of its ion 
via the hyperfine interaction. Electron spins interact with each other through the Ising 
interaction. All electron spins are in their ground states. The nuclear spins also interact with 
each other through the Ising interaction. This interaction is responsible for quantum logic 
gates. The dipole-dipole interaction is suppressed due to the use of the magic angle between 
the chain and the direction of the external magnetic field. The key point of our model 
is the following: we consider the interaction (for simplicity, the Ising interaction) between 
an electron spin and its two neighboring nuclear spins. This interaction can originate, for 
example, if the electron density at the neighboring nuclei is not zero. Thus, the electron 
spin frequency for a particular ion can take eight values. As an example, for the B-ion one 
has, 

U;''' = + Jel^ + Je^ + ^^/2 + ^el^ + (29) 

= + J^f + Jfe^ - + J^^ + 

and so on. Here cu^^'' corresponds to the electron frequency for the case in which the nuclear 
spin of the same ion is in the state the nuclear spin of the neighboring A-ion is in the 
state \k), and the nuclear spin of the neighboring C-ion is in the state ajg is the Zeeman 
frequency, J^'^ and J^J-^ are the constants of the electron-electron interaction, is the 
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hyperfine constant (in the frequency units), and J^^ and J^^ are the constants of the 
electron-nuclear interaction for neighboring ions. 

Thus, the ESR frequency depends on the position of an electron spin in the chain (because 
the Zeeman frequency is nonuniform) and the states of three nuclear spins - the nuclear spin 
of its own ion (via the hyperfine interaction) and the nuclear spins of two neighboring ions. 
One can tune the frequency of an electromagnetic pulse in such a way that it is resonant 
with the electron spin only if it is in a definite position in the chain and the three nuclear 
spins mentioned above are in the definite states. 

The strategy for phase correction for the quantum full adder is the following. The full 
adder is implemented by a sequence of the nuclear 7r-pulses. A nuclear 7r-pulse causes a 
phase shift of 7r/2 for the resonant states. There are six possible phase distortions for non- 
resonant states (3 possible states for neighboring nuclear spins times two possible states 
for the selected nuclear spin). Correspondingly, one has six possible frequencies for the 
electron spin of the selected ion in the non-resonant state. Because of the large value of the 
electron gyromagnetic ratio, we assume that each of the corresponding electron transitions 
can be driven without noticeable non-resonant effects on the electron transitions with close 
frequencies. After the nuclear 7r-pulse, one applies 12 electron 7r-pulses: two 7r-pulses for 
every possible frequency of the selected electron spin in the non-resonant state. The phase 
of the first 7r-pulse is zero, the phase of the second one is 0. The total phase shift of the 
wave function of the ion is (tt -|- 0). By choosing an appropriate value for 0, one can change 
the phase shift for a specific non-resonant state to the value of 7r/2, which corresponds 
to the resonant state. After the action of 12 electron 7r-pulses, all states in the quantum 
superposition will have the same phase, and all electron spins will be in their ground state. 
Thus, the phase distortion generated by a nuclear 7r-pulse can be corrected with 12 electron 
TT-pulses. 
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VIII. CONCLUSIONS 



We have demonstrated a strategy for implementation of a full adder in a nuclear spin 
quantum computer. To add an arbitrary superposition of L-qubit "quantum numbers" and 
a fixed L-bit "classical number" our scheme requires 2L+1 qubits and less than 27L resonant 
TT-pulses. We have simulated the action of this full adder for a spin chain containing 201 
qubits. Using the 27rfc-method allowed us to minimize the generation of unwanted non- 
resonant error states. These error states are shown to have a few preferred discrete values of 
probability (a "line spectrum"). Our simple spin model is generalized to correct undesired 
phase distortions. 
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